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This  report  comprises  two  analytical  studies.  In  the  first,  the 
methods  of  potential  theory  are  applied  to  the  problem  of  determin¬ 
ing  the  velocities  In  the  approaches  to  a  canal  connecting  an  ocean 
with  a  lagoon  or  bay.  Entrances  with  and  without  jetties  extending 
into  the  ocean  are  considered.  The  streamlines  are  determined  and 
velocities  along  the  streamlines  are  given  in  a  form  suitable  for 
simple  application  to  engineering  problems.  In  the  second  study, 
Tollmien's  analyses  for  (l)  the  free  jet  boundary  and  (2)  jet  ef¬ 
flux  from  a  linear  slot  are  applied  to  the  problem  of  the  efflux  of 
a  two-dimensional  Jet  from  a  slot  of  definite  width  with  particular 
reference  to  Jet  flow  from  a  channel  into  a  region  of  surrounding 
quiet  water,  such  as  an  ocean.  The  literature  was  searched  for  ap¬ 
propriate  experimental  data,  arid  the  constant  of  Tollmien's  theory 
is  evaluated  accordingly.  The  streamlines  are  determined,  and  ve¬ 
locities  along  the  streamlines  computed  and  presented 

in  a  manner  susceptible  of  simple  application 

to  practical  engineering  problems. 
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PREFACE 


The  Corps  of  Engineers,  Department  of  the  Army,  is  charged  hy  Con¬ 
gress  with  the  responsibility  of  establishing  and  maintaining  navigable 
channels  to  the  major  seaports  of  the  United  States.  This  responsibility 
frequently  involves  the  establishment  of  a  navigable  channel  at  tidal  in¬ 
lets  and  estuaries.  To  advise  the  Chief  of  Engineers  in  this  phase  of  his 
work,  a  Committee  on  Tidal  Hydraulics  has  been  set  up  within  the  Corps  of 
Engineers.  This  Committee,  among  other  duties,  arranges  to  have  technical 
studies  made  of  various  problems  connected  with  the  establishment  and 
maintenance  of  channels  in  tidal  waterways. 

One  of  the  problems  encountered  at  tidal  entrances  is  that  of  flow 
patterns  -  velocity  and  direction  -  which  develop  at  the  various  inlets 
selected  for  navigation  improvements.  This  problem  has  not  received  suf¬ 
ficient  attention  in  the  past  to  bring  about  a  quantitative  understanding 
of  these  flow  patterns.  As  one  step  in  overcoming  this  lack  of  under¬ 
standing,  the  Committee  requested  assistance  from  the  Mechanics  Division  of 
the  Bureau  of  Standards  in  laying  the  theoretical  groundwork  for  the  solu¬ 
tion  of  this  problem.  Arrangements  for  this  assistance  were  made  by  the 
Beach  Erosion  Board  of  the  Corps  of  Engineers  in  a  letter,  dated  20  October 
1950,  to  the  National  Bureau  of  Standards . 

The  problem,  in  effect,  has  two,  more  or  less  independent  aspects: 
the  definition  of  the  pattern  of  flow  from  the  ocean  into  the  inlet,  and 
the  definition  of  the  pattern  of  flow  from  the  inlet  into  the  ocean.  These 
two  aspects  can  be  conveniently  handled  separately,  and  the  National  Bureau 
of  Standards  reported  on  them  separately  to  the  Committee  in  1951.  The  com¬ 
mittee  has  chosen  to  publish  the  two  reports  together  in  this  volume,  but 


ill 


has  maintained  the  separation  of  the  two  studies  in  the  publication. 

The  funds  available  did  not  permit  as  extensive  an  exploration  of 
these  problems  as  is  desirable.  However,  the  two  reports  that  have  been 
prepared  are  believed  to  represent  a  valuable  contribution  to  the  litera¬ 
ture  on  this  subject.  Certain  simplifying  assumptions  were  made  in  both 
cases,  and  only  a  limited  number  of  conditions  were  investigated;  these 
are  explained  in  detail  in  the  report. 

Velocity  Pattern  at  Entrance 

The  first  part  of  the  report,  prepared  by  John  L.  French,  National 
Bureau  of  Standards,  deals  with  the  pattern  developed  by  flow  from  the 
ocean  into  an  inlet.  For  this  study,  the  following  simplifying  assumptions 
were  made: 

a.  The  depth  was  uniform  and  the  same  in  the  ocean  and  the  inlet. 

b.  The  effect  of  viscosity  was  neglected. 

Two  forms  of  inlets  were  considered:  an  inlet  without  jetties,  and  an  in¬ 
let  with  jetties  extending  into  the  ocean. 

The  author  recognized  that  the  simplifications  introduced  into  the 
analysis  would  prevent  an  exact  correlation  between  the  flow  patterns  de¬ 
veloped  by  his  study  and  actual  flow  patterns  under  field  conditions.  He 
points  out  that  this  discrepancy  would  become  more  and  more  evident  at 
points  downstream  from  the  entrance  where  diffusion  would  cause  the  jet  to 
spread  over  the  entire  width  of  the  jetty  channel. 

This  analysis  can  he  used  to  predict  the  flow  direction  and  velocity 
pattern  at  the  inlet  for  conditions  of  uniform  depth  and  unobstructed  ap¬ 
proaches.  The  presence  of  shoals  or  other  obstacles  would,  of  course, 
cause  the  actual  flow  pattern  to  vary  from  the  predicted  patterns  of 
streamlines  as  shown  in  figs.  5  and  8  of  the  report.  An  example  of  the  use 
of  this  analysis  of  flow  at  tidal  entrances  is  given  in  Appendix  A  to  this 
report . 
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Plow  Pattern  at  Exit 


The  second  part  of  the  report,  also  prepared  by  Mr.  French,  deals 
with  the  pattern  of  flow  from  an  inlet  into  the  ocean.  The  physical  as¬ 
pects  of  this  problem  suggest  similarity  to  the  diffusion  of  a  submerged 
two-dimensional  jet,  and  the  velocity  profiles  have  been  developed  on  this 
basis.  An  example  of  the  use  of  this  analysis  of  discharge  from  an  inlet 
into  the  ocean  is  given  in  Appendix  B, 
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The  author  is  indebted  to  Dr.  Garble  H.  Keulegan  of  the  national 
Bureau  of  Standards*  staff  for  the  formulation  of  the  two  problems  herein 
discussed  and  for  many  valuable  suggestions  during  the  conduct  of  the 

investigation. 
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TIDAL  FLOW  IN  ENTRANCES 


PART  I:  THE  YELOCITI  DISTRIBOTIOl  AT  THE  ENTRANCE 


Scope  of  Investigation. 

1.  At  a  conference  between  Dr.  Martin  A«  Mason  of  the  Beach  Erosion 
Board  and  Dr.  Garbis  H.  Keulegan  of  the  National  Bureau  of  Standards,  the 
scope  of  the  project  was  defined  in  detail  and  the  following  four  aspects 
of  tidal  flow  from  an  ocean  through  a  relatively  short  channel  into  a 
lagoon  were  agreed  upon  for  the  investigation: 

a.  Velocity  distribution  at  the  entrance. 

b.  Velocity  distribution  at  the  exit. 

c.  Determination  of  the  mean  velocity  in  the  channel. 

d.  Velocity  distribution  in  the  channel. 

The  Present  Problem 


2.  Part  I  of  this  report  is  confined  to  the  first  problem  listed 
above— the  velocity  distribution  of  the  water  as  it  approaches  and  enters 
the  connecting  channel  between  the  ocean  and  the  lagoon. 

3.  Fig.  1  is  a  definition  sketch  of  the  ocean,  connecting  channel, 

and  lagoon.  Inasmuch  as  the  ocean, 
connecting  channel,  and  lagoon  in 
question  are  assumed  to  be  subject 

to  tidal  action,  this  aspect  of  the  - - - — 

investigation  will  apply  to  flow 

from  the  lagoon  to  the  ocean  as  lagoon 

well  an  flow  from  the  ocean  to  the  |  channel 

lagoon.  The  various  factors  af¬ 
fecting  the  velocity  distribution 
at  and  near  the  entrance  include, 
in  addition  to  the  shape  or  geom¬ 
etry  of  the  entrance  itself,  the 
variation  in  depth  of  the  ocean  or 


Fig.  1.  Definition  sketch  of  problem 


template!  as  forming  a  part-  of  the  project.  For  the  present,  the  effect 
of  frictional  resistance  and  the  effect  of  changing  depth  in  the  ocean  or 
lagoon  as  the  channel  entrance  is  approached  will  he  neglected.  Hence  the 
present  problem  resolves  itself  into  the  effect  of  various  channel-entrance 
shapes  on  the  velocity  distribution  of  a  two-dimensional  flow  of  a  perfect 
fluid  from  a  large  body  of  water  into  a  relatively  narrow  channel  entrance. 

k.  The  velocity  picture  at  the  entrance  will  be  determined  by  the 
methods  applicable  to  free  streamlines,  methods  devised  by  Helmholtz  and 
Kirchhoff  and  later  elaborated  by  various  investigators.  If  fluid  parti¬ 
cles  in  flowing  around  a  sharp  corner  of  radius  zero  are  assumed  to  remain 
in  contact  with  the  bounding  surface,  an  infinite  acceleration  of  the  fluid 
particles  is  required  to  accomplish  this  result.  This  in  turn  requires 
that  the  velocity  and  pressure  he  infinite  at  the  corner.  In  order  to 
overcome  this  difficulty,  Helmholtz  introduced  the  concept  of  free  stream¬ 
lines,  in  which  the  assumption  is  made  that  separation  occurs  in  the  flow 
of  fluid  particles  around  sharp  comers.  Thus,  downstream  from  abrupt 
changes  in  the  bounding  surface,  such  as  those  considered  in  this  report, 
the  streamlines  of  the  flow  are  separated  from  the  boundary  by  an  area  oc¬ 
cupied  by  fluid  assumed  to  be  at  rest. 

5.  With  real  fluids,  owing  to  viscosity,  the  motion  is  so  modified, 
of  course,  that  points  of  infinite  velocity  and  pressure  do  not  occur,  nor 
does  the  fluid  between  the  free  streamline  and  the  bounding  surface  remain 
at  rest;  it  is  in  eddying  motion.  With  the  assumption  of  free  streamlines, 
the  velocity  and  hence  the  pressure  at  abrupt  corners  is  finite,  and  to 
this  extent  the  solutions  by  this  method  of  various  flow  problems  involving 


of  separation.  The  region  between  the  free  streamline  ana  the  channel 
wall  is  in  eddying  motion,  and  the  lateral  diffusion  of  the  eddies  as  they 
move  downstream  obliterates  the  potential  cor®  of  the  jet  and  causes  the 
outside  jet  boundaries  to  expand  to  the  full  -width  of  the  channel  at  a 
distance  on  the  order  of  10  channel  widths  from  the  entrance. 

6.  Under  these  circumstances,  the  use  of  the  free  streamline  method 
to  determine  velocity  distributions  at  distances  in  excess  of  approximately 
two  or  three  channel  widths  downstream  from  the  entrance  would  not  be 
warranted. 

7.  The  solutions  in  this  report  are  given  in  graphical  fora,  i.e., 
the  streamlines  are  shown  and  the  velocities  along  the  streamlines  are 
given  in  terms  of  the  free  streamline  velocity  and  also  in  terms  of  the 
mean  velocity  in  the  channel.  The  mean  velocity  in  the  channel  must  be  de¬ 
termined  from  a  consideration  of  the  tidal  fluctuations  in  the  ocean  and 
lagoon,  and.  the  dimensions  and  roughness  of  the  connecting  channel.  This 
problem  will  be  considered  in  a  later  report  on  this  project. 


Previous  Treatment  of  Free  Streamline  Flow 


The  literature  is  replete  with  investigations  of  two-dimensional 
bleats  in  which  the  bounding  surfaces  take  the  form  of  possible 
entrance  shapes.  Unfortunately ,  however,  these  problems  have  been 
xeluslvely  considered  from  the  standpoint  of  determining  the  equa- 
the  free  streamlines  and  the  coefficient  of  contraction.  The  com 
of  the  equipotential  lines  and  streamlines  in  the  body  of  the 
fluid  has  been,  in  the  great  Majority  of  cases,  entirely  omitted. 


nave  application  to  certain  tidal-entrance  shapes  is  s 
the  references  cannot  be  individually  described  herein 
vestigatlons  have  been  so  numerous  that  although  an  ex 
literature  was  made,  it  was  by  no  means  exhaustive,  an 
in  view  of  the  time  limitation  on  the  project.  Nevertl 
resume  of  the  more  important  references  is  believed  wa: 

10.  First,  most  of  the  standard  references  in 


treat  the  free  streamline, 

through  an  aperture  in  a  plane  wall,  and  through  a  channel  corresponding 

to  the  two-dimensional  form  of  Borda*s  mouthpiece.  Without  exception,  tl 

treatment  is  from  the  viewpoint  of  determining  the  equation  of  the  free 

streamline  and  determining  the  coefficient  of  contraction  of  the  jet. 

7 

11.  Prasil  gives  a  solution  of  the  two-dimensional, 

flow  through  an  aperture  in  a  plane  wall  and  shows  a  flow  net. 

the  scale  of  the  coordinate  system  is  not  given,  and  other  inadequacies 

limit  the  usefulness  of  the  flow  net  given. 

2 

12.  Greenhill  has  solved  a  variety  of  problems  in  two-dimensional 
flow  which  have  application  to  possible  channel  entrances  of  various  shapes. 
Among  these  are  those  shown  in  fig.  2.  The  angle  0  in  fig.  2(a)  and 

2(h)  may  he  varied  to  give  solu¬ 
tions  to  many  special  problems, 

among  them  being  most  of  those 

9 

treated  by  von  Mises.' 

5 

13.  Michell  also,  among 
many  others,  solved  the  problem  of 
the  modified  Borda’s  mouthpiece 
shown  in  fig.  2(d). 

O 

lh.  Greenhilv  using  pri¬ 
marily  the  work  of  Cisotti,  Levi- 
Civita,  Levy,  and  others,  has 
given  solutions  to  many  two- 
dimensional  flow  patterns  involv¬ 
ing  curved  boundaries, 
these  are  shown  in  fig. 


M 


W 


(d) 

Two-dimensional  flow 
boundaries 


15.  Cisotti  has  also  published  solutions  to  most  of  the  flow 
problems  shown  in  figs.  2  and  3. 

16.  The  references  cited  above  give  solutions  to  the  various 
problems  indicated,  generally  by  giving  the  functional  relation  between 
the  complex  potential  and  a  parameter,  which  in  turn  is  related  through 


*  Raised  numbers  refer  to  references  at  the  end  of  this  part  (Part  i). 
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another  factional  expression  to  the 
physical  z  plane.  The  actual  de- 


17.  In  general,  the  work  in- 

(<A  (*>} 

volved  in  applying  the  foregoing 
solutions  to  the  determination  of  the 
flow  net  is  not  difficult!  neverthe¬ 
less,  it  is  in  most  cases  tedious,  \ 

and  the  amount  of  labor  involved  is 

by  no  means  small.  For  this  reason,  \  / 

it  has  been  impossible  with  the  tine 

(c] 

at  our  disposal  to  utilise  all  the 

material  available  and  to  compute  Fig.  3*  Two-dimensional  flow 

the  flow  net  for  all  the  boundary  wtth  curved  boundaries 

forms  shown  in  figs.  2  and  3*  I® 

the  remainder  of  this  report,  the  equipotential  and  streamlines  for  the 
two-dimensional  flow  through  an  aperture  in  a  plane  wall  and  through  the 
two-dimensional  form  of  a  Borda's  mouthpiece  have  been  determined.  Inas¬ 
much  as  any  streamline  may  he  replaced  by  a  hounding  surface  without  in 
any  way  disturbing  the  flow  pattern,  these  two  examples  will  yield  a  fairly 
representative  group  of  possible  channel  entrances.  It  is  nevertheless 
regretted  that  more  of  the  examples  cited  could  not  have  been  applied  to 
the  present  problem,  and  it  is  believed  that  the  future  application  of 
these  two-dimensional  flow  solutions  to  the  determination  of  the  equi¬ 
potential  and  streamlines  for  the  various  shapes  of  possible  channel  en¬ 
trances  will  prove  a  useful  tool  in  the  practical  solution  of  the  tidal- 


termination  of  ¥eloc:i 


water..  Since  this  difference  in  elevation  is  caused  by  tidal  action,  the 
flow  Into  the  channel  will  vary  with  time.  However,  since  the  tidal  action 
will  vary  slowly  with  time,  we  can  for  all  practical  purposes  consider 
the  problem  of  the  velocity  distribution  at  the  channel  entrance  as  being 
one  involving  steady  flow  only. 

19.  Owing  to  the  fact  that  the  velocity  along  a  streamline  will  vary 
as  the  channel  entrance  is  approached,  it  is  obvious  from  a  consideration 
of  Bernoulli  s  equation  that  the  surface  elevation  will  vary.  However,  the 
velocities  that  are  encountered  in  tidal  channels  are  such  that  the  varia¬ 
tion  in  surface  elevation  due  to  variation  in  velocities  will  he  small  com¬ 
pared  to  the  depth  of  the  channel.  Therefore,  insofar  as  surface  varia¬ 
tions  are  concerned,  the  problem  may  he  considered  for  practical  purposes 
to  he  essentially  two-dimensional.  As  noted  previously,  the  effects  of 
friction  and  variation  in  depth  of  the  ocean  as  the  channel  entrance  is  ap¬ 
proached  will  he  treated  In  a  subsequent  report,  and  for  the  present,  con¬ 
stant  depth  and  nonviscous  flow  are  assumed. 

20.  Under  these  conditions,  the  methods  of  conformal  mapping  have 
direct  application  to  the  problem  of  obtaining  the  velocity  distribution 
at  the  channel  entrance.  For  irrotatlonal,  two-dimensional,  steady-flow 
problems  it  is  necessary  only  that  laplace's  equation  in  two  dimensions  he 
satisfied,  and  that  the  boundary  conditions  regarding  coincidence  of 
streamlines  with  fixed  boundaries  be  observed.  Where  there  is  a  free 
streamline,  the  further  requirement  Is  made  that  the  pressure  and  hence 
the  velocity  along  the  free  streamline  he  constant. 

21.  The  ability  of  the  processes  of  conformal  transformations  to 
satisfy  the  above-mentioned  requirements  and  hence  to  provide  solutions  to 
irrotatlonal,  two-dimensional  flow  problems  has  been  amply  demonstrated  in 
most  of  the  standard  references  on  hydrodynamics,  and  there  is  no  need  to 
repeat  the  argument  herein. 

Entrance  analogous  to 
two-dimensional  orifice 

22.  This  is  the  type  of  entrance  from  the  ocean  to  the  channel  (and 
from  the  lagoon  to  the  channel)  shown  in  fig.  1. 
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Greentd.il  (page  33  of  reference  2}  bas 
shown  for  the  flow  orientation  given  in 
fig.  4  that  the  relation  between  the 
flow  pattern  in  the  physical  plane 
z  =  x  +  i  y  of  the  flow,  and  the  com¬ 
plex  potential  ¥  =  $  +  if,  where  0 
and  i|r  are  respectively  the  potential 
and  stream  functions,  is  given  by 

1/2  f  -  log  - 1  u> 


and 


sinh  Q  a  e  (2) 

where  £  is  related  to  ¥  ,  0,  ,  and 
other  characteristics  of  the  flow 
through  the  relation 


y 


Pig.  4.  Channel  entrance 
corresponding  to  a  two- 
dimensional  orifice 


dz  _  v  &±e  _  n 

dW 


(3) 


Tn  the  preceding  expressions,  c  is  the  half -width  of  the  jet  at  ini  ini  fcy 
as  shown  in  fig.  4;  V  is  the  velocity  along  the  free  streamline  and  also 
the  velocity  of  the  jet  at  infinity;  9  is  the  angle  that  the  velocity 
vector  makes  with  the  positive  x-axis;  Q  is  the  rate  of  discharge  per 
unit  depth  of  channel  and  therefore  has  the  dimensions  of  a  velocity  multi¬ 
plied  by  a  length;  and  a  Is  the  velocity  at  any  point  x  ,  y  in  the  flow 
pattern.  The  potential  and  stream  functions  j>  and  *  are  related  to  the 
x  and  y  velocity  components  through  the  expressions 


and 


^  ^  ^  ¥  SSS 


velocity  multiplied  by  a  length  so  that  j#?/Q  and  consequently  £  and  ft 
are  dimensionless. 

24.  Equation  1  has  been  derived  on  the  basis  that  the  zero  potential 

line  passes  through  the  two  points  B  and  B*. 

25.  In  equation  2,  replacing  ft  by  its  value  ft  =  log  £  from 
equation  3,  and  replacing  the  hyperbolic  sine  by  its  exponential  form, 
we  have 


which  yields 


(4) 


and  since 
sign  and 


5  is  infinite  when  W  is  infinite,  we  must  use  the  positive 


rij  f~2 iflf 

t;  =  eQ  +  le  Q  +  1 


(5) 


Or  we  may  write 


i  = 


I  ~~2rfW 

e  Q  ll  +  f  1  +  e”  Q 


(6) 


26.  The  real  and  imaginary  parts  of  £  may  he  conveniently  deter¬ 
mined  by  replacing  the  terms  under  the  radicals  in  equations  5  and  6  by 
their  equivalent  series  forms.  The  binomial  expansion  equivalent  to  the 
radical  term  in  equation  5  will  he  convergent  for  negative  values  of  the 
potential  0  ,  and  equation  6  may  he  used  to  compute  £  for  positive 
values  of  0  •  Therefore,  letting  I  and  tj  he  the  real  and  imaginary 
parts  of  £  ,  determined,  from  equations  5  we  have 


(7) 


9 


(8) 


(9) 


28.  letting  b  be  the  half -width  of  the  channel  as  indicated  in 
fig.  k,  we  have  y  «  -b  at  $  -  0  and  if  *  0/2*  from  which  by  equation  4, 

5  .  i  at  the  point  B'.  Substitution  in  equation  9  gives 


£  _  -  (10) 

b  ~  2  +  « 


29.  Equations  7,  11,  and  12  permit  the  computation  of  the  coordi¬ 
nates  in  the  z  plane  corresponding  to  any  value  of  the  complex  potential 
¥  .  These  equations  have  been  used  to  compute  values  of  x/b  and  y/b  cor¬ 
responding  to  various  values  of  *0/Q  and  t/Q.  The  results  are  given  in 
table  1.  In  fig.  5  the  data  of  table  1  have  been  used  to  plot  the  stream¬ 
lines  and  equipotential  lines  of  the  flow  net.  Comparison  of  fig.  5  with 
fig.  4  will  show  that  the  orientation  differs  by  l80°.  This  was  done 

merely  for  convenience  in  drawing  fig.  5- 

30.  Since  the  velocity  is  tangential  to  the  streamlines  at  all 
points,  the  data  of  fig.  5  indicate  the  direction  of  the  current  as  it  ap¬ 
proaches  the  channel  entrance.  The  problem  of  determining  the  magnitude 
of  the  velocity  remains. 

Ti. ,  Since  a  *  8^/cls,  where  q  is  the  velocity  measured  in  the 


(15) 


Eliminating  0  between  equations  l4  and,  15  and  rearranging. 


F  2s£  \2 

2it0 

cosh*  log  ~  =  l/2 
<1 

0  ^  Hh  X  Hb  ^ 

-  Ite  «  sin2  d 
<4 

Since  v/Q  must  be  infinite  for  0  =  <*>,  the  positive  sign  before  the  radical 
must  be  used. 

32.  Through  this  equation  the  velocity  of  any  point  may  be  deter¬ 
mined  in  terms  of  the  constant  velocity  ¥  of  the  free  streamline.  In 
table  2  this  has  been  done  for  various  values  of  t/Q  and  «0/Q.  It  is, 
however,  more  convenient  in  many  practical  applications  to  express  q  in 
terms  of  the  mean  velocity  in  the  channel.  Letting  V  be  the  mean  channel 
velocity,  we  have 


where,  as  before,  c  is  the  half -width  of  the  jet  at  infinity  and  b  is 
the  half -width  of  the  channel.  In  view  of  equation  10, 

y  ,  y  «  0.611  V 

m  2  +  it 

33.  Using  this  relation  and  equation  l6,  values  of  q/V  and  q/v 
have  been  plotted  against  s/b  in  fig.  6  where  s  is  the  distance  measured 
along  a  particular  streamline  from  its  intersection  with  the  zero  po¬ 
tential  line.  The  value  of  s  for  any  given  value  of  0  and  f  was  de¬ 
termined  from  fig,  5* 

3k .  The  use  of  the  data  of  figs,  5  and  6  in  determining  the  hori¬ 
zontal  velocity  distribution  in  the  approaches  to  a  channel  entrance  are 
illustrated  by  the  following  example.  Suppose  the  half -width,  b  ,  of  tbe 
channel  is  500  ft  and  the  difference  in  surface  elevation  between  the 
ocean  and  lagoon  is  such  that  the  rate  of  discharge  in  the  channel,  of 
width  2b,  is  2000  cfs  per  foot  of  depth.  Suppose  further  that  we  wish  to 
determine  the  velocity  at  the  point  x/b  =  0,  y/b  =  0.7,  or  in  this  example 
for  a  channel  half -width,  b  ,  of  500  ft,  x  »  0  and  y  =  350  ft.  From 


Fig.  6.  Variation  of  velocity  along  streamlines  for  flow  through 

two-dimensional  orifice 

fig.  5  this  point  is  seen  to  he  on  the  streamline  for  which  f/<£  =  3/8,  an 
using  fig.  5,  s  measured  from  the  zero  potential  line  to  the  point  in 
question  is  found  to  be  -0.3b.  Referring  to  the  curve  for  ty/Q  =  3/8  in 
fig.  6,  it  is  found  that  q/Vffi  =1.26  for  the  s/b  =  -0.3.  Since  V  =  Q/2b 
we  have,  using  the  value  of  Q  and  b  assigned  for  this  illustration, 

V  =2.00  ft/sec 

Hence  at  s/b  =  -0.3,  on  the  streamline  in  question  we  have 

q  =  1.26  Vffl  =  2.52  ft/sec 

35.  The  direction  of  the  velocity  vector  will,  of  course,  he 
tangential  to  the  streamline  at  the  point  in  question. 

36.  An  example  of  the  application  of  the  generalized  solution  show 
in  fig.  5  to  a  specific  assumed,  condition  is  given  in  Appendix  A. 

3  7 *  It  is  again  pointed  out  that  the  analysis  from  which  the  pre¬ 
ceding  results  are  derived  is  based  on  the  assumption  that  pure,  nonvisco 
potential  flow  obtains,  and  that  the  jet  shown  in  fig.  k  is  bounded  by  fr 
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walls  being  filled  with  water  at  rest.  In  reality  the  region  between  the 
free  streamline  and  the  channel  wall  is  in  eddying  motion,  and  the  lateral 
diffusion  of  the  eddies  as  they  move  downstream  causes  the  jet  to  spread 
eventually  to  the  full  width  of  the  channel.  Under  these  circumstances,  as 
previously  stated,  the  use  of  the  results  herein  obtained  to  determine  ve¬ 
locity  distributions  at  appreciable  distances  downstream  from  the  channel 
entrance  would  not  be  justified. 

38.  Since  the  velocity  component  normal  to  a  streamline  is  zero,  any 

of  the  streamlines  of  fig.  5  may  be  replaced  by  a  bounding  surface  without 

in  any  way  affecting  the  flow  net.  In  this  manner  a  relatively  wide 

variety  of  possible  channel-entrance  conditions  may  be  obtained  from  fig.  5* 

In  this  connection  it  is  obvious  that  the  relative  magnitudes  of  the  mean 

velocity  in  the  channel,  and  ¥  the  velocity  along  the  free  streamline, 

and  also  the  velocity  of  the  jet  at  relatively  large  distances  from  the 

channel  entrance  depend  on  the  type  of  channel  entrance  assumed.  For 

example,  with  the  channel  width  being  assumed  equal  to  the  full  width  of 

the  two-dimensional  orifice,  ¥  **  0.6ll  ¥  as  we  have  seen.  However,  if  we 

m 

replace  the  free  streamline  with  a  fixed  boundary,  c  =  b  and  =  ¥. 
Obviously  the  values  of  q/V  given  in  fig.  6  apply  only  to  the  case  of  the 
channel  shown  in  fig.  5» 

39.  Since  the  free  streamlines  approach  asymptotically  the  line 
y/b  =  0.611,  the  channel  formed  by  replacing  the  free  streamlines  by  a 
fixed  boundary  will  not,  strictly  speaking,  consist  of  a  channel  with 
straight  parallel  walls.  However,  it  is  apparent  from  fig.  5  that  the 
free  streamlines  approach  the  value  y/b  =  0.6ll  so  rapidly  for  small  values 
of  x/b  that  a  channel  consisting  of  straight  parallel  walls  beyond 

x/b  =  -1,  and  with  the  walls  coinciding  with  the  free  streamlines  from 
x/b  »  -1  to  x/b  =  0,  may  for  all  practical  purposes  be  considered  to 
follow  the  free  streamlines  throughout  their  length. 

40.  The  modification  in  the  method  of  computing  the  velocity  at  any 
point  by  replacing  a  streamline  of  fig.  5  hy  a  fixed  boundary  is  obvious. 
For  example,  to  continue  the  illustration  of  replacing  the  free  streamline 
with  a  bounding  surface,  we  see  from  fig.  6  that  q /v  =  0.77  at  the  point 
on  the  streamline  Y/O.  —  3/8  £«r  which  s/b  =  -0.3* 


i4 


4.1 .  Since  ¥  =  V_t  =  Q/2b  in  this  case,  using  the  values  of  the 
previous  illustration,  we  have 

¥  =  1000/500  =2.0  ft/sec 


and 


q  -  0.T7  ¥  =  1.5^  ft/sec 

42.  If  hoth  boundaries  of  the  flow  are  taken  to  be  streamlines  in 

the  interior  of  the  jet,  it  is  again  evident  from  fig.  5  that  a  channel 

with  parallel  straight  walls  for  x/b  <  -1  and  following  the  streamlines  in 

question  for  s/b  >  -1  closely  approximates  the  theoretical  configuration, 

and  that  again  as  for  the  case  of  the  free  streamline  being  replaced  by  a 

fixed  wall,  c  is  to  a  close  approximation  equal  to  b  ,  and  likewise  from 

a  consideration  of  fig.  6  it  is  evident  that  ¥  may  with  small  error  he 

m 

taken  as  ¥  .  Hence,  values  of  q/v  shown  in  fig.  6  may  be  used  to  obtain 
the  velocity  at  any  point  in  the  field  of  flow  when  the  channel  and  en¬ 
trance  are  such  that  the  boundary  walls  consist  of  straight  parallel  walls 
up  to  x/b  =  -1  and  coincide  with  streamlines  as  x/b  Increases  in  the  posi¬ 
tive  direction. 

43.  It  will  be  noted  in  fig.  5  that  the  stream  function  if  has  been 
expressed  in  terms  of  Q  ,  the  discharge  per  unit  depth  of  channel.  This 
appears  logical  since  it  is  a  characteristic  of  the  stream  function  that 
the  difference  between  its  value  at  any  two  points  represents  the  volume 
rate  of  flow  per  unit  depth  across  any  line  joining  the  two  points.  How¬ 
ever,  the  stream  function  may  he  expressed  equally  well  in  other  forms. 

For  example,  since  Q  =  2  ¥c,  and  using  equation  10,  the  lines  of  constant 

ijf  may  he  conveniently  represented  by  t/vb  »  7,  where  7  is  a  dimension¬ 
less  number. 

Long  parallel  jetties  extending 
into  ocean  from  the  channel  entrance 

44.  This  type  of  entrance  condition  corresponds  essentially  to  the 
two-dimensional  form  of  Horde. ’ s  mouthpiece — the  longer  the  jetties,  the 
closer  the  correspondence. 


45.  In  fig.  7  is  shown  the  geom¬ 
etry  of  the  channel  entrance.  The 
subscript  <®  on  some  of  the  points  is 
used  to  indicate  that  the  point  is 
located  at  infinity.  As  before,  Q 
is  the  discharge  into  the  channel  per 
unit  depth  of  channel,  b  is  the 
half -width  of  the  channel,  and  the 
streamlines  will  he  denoted  in  terms 
of  the  discharge  per  unit  depth  of 
channel.  Thus,  the  streamline  along 
the  boundary  ,  A*,  c;  being  *  =  -Q/2  and  the  upper  bounding  streamline 

being  t  «  q/2,  the  zero  streamline  is  coincident  with  the  x-axis. 

46.  Milne-Thompson  (page  283  of  reference  6),  proceeding  by  the 
usual  methods  of  conformal  transformations,  derives  the  solution 


-v§4  +  V?T"i]2 

(17) 

¥  _  2  °  bV  log  t  -  i  a  hV 

(18) 

where  t  is  a  complex  parameter;  b  is  the  half -width  of  the  channel  as 
shown  in  fig.  7;  0  is  the  coefficient  of  contraction;  and  ¥  is,  as  be¬ 
fore,  the  velocity  along  the  free  streamline.  Obviously  Q  ,  the  discharge 
per  unit  depth  in  the  channel,  is  equivalent  to  the  term  2  a  bV  and 

equation  IS  becomes 

W  =  ~  log  t  -  i  ■§ 

It  et 

(19) 

from  which 

Wtf  A  lit  Wn 
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cu 
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Fig.  7.  Channel  entrance  corre¬ 
sponding  to  long  parallel  jetties 
extending  into  ocean 


Substituting  the  above  relation  in  equation  17 
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2iM 

e  Q  +  1  +  S  log 

it 


*W 


+  i 


+  H  4-  A 


where  A  may  be  a  complex  constant*  In  order  to  evaluate  A  we  faave^ 
since  in  the  development  of  equation  IB  the  potential  0  was  assumed  to  "be 
zero  at  A  and  A", 


at  z  =  lb,  W  =  i.Q/2 


at  2  =  -lb,  W  ■  -i.Q/2 


The  application  of  the  first  condition  to  equation  22  yields 


A  =  IVb  -  iQ  +  Q/« 


(23) 


and  the  second  condition  gives 


A  «  -iVb  +  iQ  +  Q/  it 


(24) 


Hence  A  -  Q/it  and  equation  22  becomes 


2*W 


Vz  =  ~  e  Q  +  —  e 

it  it 


m 
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+  1  +  “7  log 
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sW  /  2kW 
e  + 


+  1 
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From  equations  23  and  24  it  is  clear  also  that 


Vb  m  Q  (26) 

and  since  Q  =  2cV,  where  c  is  the  half -width  of  the  jet  at  infinity,  it 
follows  that  the  coefficient  of  contraction 

0  =  c/b  «*  l/2 


in  accordance  with  the  usual  result. 


this  connection  the  terms  under  the  radical  signs  offer  the  only  diffi¬ 
culty,  These  terms  may  he  expanded  in  a  binomial  series  as  in  the  prec 
ing  section,  or  they  may  he  evaluated  directly  in  the  usual  manner  as 
follows.  Consider  first  the  term 


appearing  in  the  right-hand  side  of  equation  25*  We  may  write 
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Placing  the  right-hand  member  of  this  relation  in  polar  form,  and  rear 
ranging,  we  have  for  positive  values  of  the  stream  function 
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-we  have  from,  equation  29 
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For  the  expression  in  brackets  appearing  in  equation  25  we  may  write 
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50.  The  two  relations  above  may  be  used  to  determine  the  equipo- 
tential  lines  and  streamlines  of  the  flow  pattern  in  the  z  plane.  In 
table  3,  values  of  x/b  and  y/b  corresponding  to  various  values  of  t/Q  and 
«0/ Q  have  been  computed;  and  in  fig.  8  these  data  have  been  used  to  plot 
the  flow  net  for  the  upper  half  of  the  channel.  Since  the  flow  is  sym¬ 
metrical  about  the  channel  axis,  the  flow  net  in  the  lower  half  will  be 
identical. 

51.  Th.e  determination  of  the  velocities  at  any  point  in  the  channel 
or  in  the  approaches  to  the  channel  may  be  made  in  a  manner  similar  to  that 
employed  in  the  preceding  section.  From  equation  21,  we  have  in  view  of 
equation  2, 
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and  using  equation  32  we  nay  write 
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Fig.  9.  Variation  of  velocity  along  streamlines  for  the  case  of  long 
parallel  jetties  extending  into  the  ocean  from  the  channel  entrance 

before  s  is  the  distance  along  the  streamline  measured  from  its  intersec¬ 
tion  with  the  zero  potential  line.  In  table  4,  computed  values  of  q/v  and 
measured  values  of  s/b  obtained  from  fig.  8  are  given.  The  data  of  figs.  8 
and  9  make  it  possible  to  compute  the  magnitude  and  direction  of  the  veloc¬ 
ity  at  any  point  in  the  field  of  flow. 

53.  It  may  again  be  noted  that  in  the  analysis  leading  to  equations 
17  and  3B,  pure  potential  flow  of  a  nonviscous  fluid  was  assumed.  In 
particular,  as  previously  pointed  out,  the  theory  assumes  that  an  abrupt 
discontinuity  in  velocity  occurs  at  the  free  streamline  and  that  the  area 
between  the  jet  and  the  boundary  is  occupied  by  liquid  at  rest.  Owing  to 
viscosity,  this  assumption  obviously  introduces  an  approximation,  and  the 
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replaced  by  a  fixed  boundary  without  in  any  way  affecting  the  flow  pattern. 
For  this  reason  the  velocity  distribution  for  a  relatively  wide  variety  of 
possible  channel-entrance  shapes  may  he  obtained  from  figs.  8  and  9* 

Status  of  Problem 


55.  In  this  part  (i),  the  flow  patterns  for  two  shapes  of  channel 
entrances  have  been  computed.  Inasmuch  as  any  of  the  streamlines  may  he 
replaced  by  fixed  boundaries  without  affecting  the  flow  pattern,  it  is  pos¬ 
sible  from  the  data  given  to  approximate  several  other  shapes  of  channel 
entrances.  It  is  nevertheless  regretted  that  the  flow  net  for  other  pos¬ 
sible  types  of  channel  entrances  could  not  he  determined  because  of  time 
limitations.  This  is  particularly  true  with  respect  to  the  channel- 
entrance  shapes  shown  in  figs,  2(a)  and  2(d)  and  also  with  respect  to  a 
solution  given  by  Michell5  which  corresponds  to  the  effect  of  a  littoral 
current  on  the  flow  through  an  aperture  in  a  plane  wall. 

56.  The  ability  of  the  data  given  to  determine  approximately  the 
velocity  distribution  in  the  approaches  to  an  actual  channel  entrance  has 

not  been  investigated.  In  this  connection  the  U.  S.  Coast  and  Geodetic 

ft 

Survey  has  in  recent  years  published  material  on  the  tidal  currents  in 
many  of  the  principal  harbors  of  the  country.  This  material  has  not  yet 
been  examined  in  detail.  If  examination  of  the  Coast  and  Geodetic  Survey 
data  should  indicate  its  suitability  for  this  purpose,  one  of  the  sub¬ 
sequent  progress  reports  of  this  project  will  apply  the  material  given 
herein  to  an  appropriate  channel  entrance  covered  by  the  Coast  and  Geodetic 
Survey  investigations. 

57.  In  this  same  connection  it  may  be  observed  that  the  boundary 
contours  of  a  natural  channel  entrance  will  in  general  never  conform  pre¬ 
cisely  to  the  simple  shapes  assumed  herein  or  indeed  to  any  other 
theoretical  investigation.  The  degree  of  approximation  introduced  thereby 
will  ordinarily  he,  at  best,  of  the  same  order  of  magnitude  as  that  intro¬ 
duced  by  the  application  of  uniform-flow  formulas  to  the  flow  of  natural 
rivers.  Under  these  circumstances,  the  degree  of  refinement  warranted  in 
theoretical  considerations  of  the  problem  is  open  to  question,  neverthe¬ 
less,  it  is  believed  that  the  extension  of  the  present  investigation  to 
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the  determination  of  the  velocity  distributions  in  some  of  the  more  com¬ 
plicated  channel-entrance  shapes,  such  as  those  of  figs.  2  and  3,  would  be 
fully  warranted. 
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PROM  A  CHANNEL  INTO  AH  OCEAN  OR  MGOON 

Introduction 

58.  In  Part  I  of  this  report,  the  velocity  distribution  in  the  ap¬ 
proaches  to  a  channel  entrance  was  considered  for  the  case  of  flow  from  the 
ocean  or  lagoon  into  the  channel  entrance.  In  Part  II,  the  reverse  problem 
will  be  considered.  That  is,  the  velocity  distribution  in  the  approaches 
to  the  channel  entrance  will  be  investigated  for  the  case  of  flow  from  the 
channel  into  the  quiet  water  of  the  ocean  or  lagoon. 

59.  A  schematic  sketch  of  the  problem  is  shown  in  fig.  10.  In  the 
figure  as  drawn,  it  is  assumed 
that  the  tidal  fluctuation  in 
the  ocean  at  the  time  of  obser¬ 
vation  is  such  that  flow  is 
from  the  lagoon  into  the  chan¬ 
nel  and  thence  to  the  ocean. 

Owing  to  turbulent  mixing  of 
the  high-velocity  channel  flow 
with  the  surrounding  quiet  Fig.  10.  Definition  sketch  of  problem 

water,  the  jet  issuing  from  the 

channel  into  the  ocean  spreads  laterally,  with  a  consequent  decrease  xn 
velocity.  If  the  ocean  depth  at  the  channel  entrance  is  greater  than  the 
depth  of  the  channel,  or  if  the  ocean  depth  increases  with  distance  away 
from  the  channel  entrance,  the  processes  of  turbulent  diffusion  will  cause 
the  jet  to  spread  vertically  from  the  channel  entrance  toward  the  ocean 
bottom.  The  complexity  of  this  three-dimensional  problem  is  such  that  it 
is  not  amenable  to  simple  analysis,  ana  for  convenience  the  present  in¬ 
vestigation  win  be  restricted  to  those  situations  in  which  the  bottom  of 
the  channel  and  ocean  are  such  that  to  an  approximation,  at  least,  the  flow 


2  6 


Statement  of  Problem 


£  IQ* 

0,  Tollmien  was  the  first  to  investigate  theoretically  the 
turbulent  diffusion  of  two-dimensional  jets.  Tollmien* s  analysis  was 
based  on  the  momentum-transport  theory  of  turbulence.  later  investigators 
have  (a)  modified  To linden's  momentum- transfer  methods  somewhat  to  secure 
better  conformance  with  experiment8;  (b)  used  the  method  of  vorticity 
transport  to  determine  the  velocity  distribution  in  the  jet6;  or  (c)  have 
proceeded  entirely  on  the  basis  of  general  momentum  and  energy  considera¬ 
tions,  with  the  use  of  the  error  function,  to  represent  the  velocity 
distribution. ” 


61.  With  suitable  choice  of  constants,  both  the  momentum-  and 
vorticity- transfer  theories  yield  velocity  distributions  in  close  agree¬ 
ment  with  experiment . 

62.  The  velocity  in  the  jet  is  also  closely  approximated  by  the 
error  function;  the  use  of  the  error  function  for  this  purpose  has  the  ad¬ 
vantage  of  simplicity,  hut  it  is  not  based  on  any  particular  theory  of 
turbulence  diffusion. 

63.  The  velocity  distributions  of  'tollmen's  jet  analysis  agree 
closely  with  experiment,  and  in  lieu  of  a  more  satisfactory  theory  of  jet 
diffusion,  the  results  of  Tollmien's  investigation  will  be  used  for  the 
present  problem  of  the  turbulent  diffusion  of  two-dimensional  jets,  in 
this  connection,  Tollmien  studied  two  aspects  of  the  problem.  The  first 
was  the  free  jet  boundary,  shown  schematically  in  fig.  11,  in  which  the 

turbulent  mixing  along  the  boundary  between 
a  stream  of  uniform  velocity  and  adjacent 
still  fluid  is  investigated.  The  second 
aspect  of  the  two-dimensional  jet  problem 
studied  by  Tollmien  was  the  case  of  jet  ef¬ 
flux  from  a  linear  slot  in  a  plane  wall,  such 
as  shown  in  fig.  12. 

64.  These  two  concepts  of  the  physical 
phenomena  involved  in  the  turbulent  diffusion 
^ jets  may  be  applied  readily  to  the 
*  Raised  numbers  reTerTxTreF^^  part  (fgrF  if)" - 


y 


Pig.  11.  Definition  sketch 
of  free  jet  boundary 


Fig,  12.  Definition  sketch  of 
TOll-mien's  jet  flow  from  a 
linear  slot 


with  uniform  velocity,  the  velocity 
discontinuity  "between  the  jet  and 
the  layer  of  still  fluid  on  either 
side  of  the  jet  generates  eddies, 
which  causes  a  diffusion  pattern  to  he 

set  up  along  the  jet  boundary  similar  to  that  shown  in  fig.  11.  The  lat¬ 
eral  mixing  process  spreads  inward  into  the  constant-velocity  core  of  the 
jet  and  outward  into  the  surrounding  still  fluid.  Consequently,  the  width 
of  the  constant-velocity  core  of  the  jet  decreases  with  distance  from  the 
slot  as  shown  in  fig.  13,  and  eventually  becomes  zero  as  the  mixing  process 

diffuses  from  both  boundaries  of  the 


two-dimensional  jet  to  the  center  line. 
This  region  of  flow  in  which  the 
constant-velocity  core  of  the  jet  is 
annihilated  has  been  aptly  termed  the 
zone  of  flow  establishment  by  Albertson, 
Dai,  Jensen,  and  Rouse This  zone 
corresponds  to  the  first  aspect  of  the 
problem  of  free  turbulence  treated  by 
Tollmien — that  is,  the  turbulent  mixing 
along  the  boundary  between  a  stream  of 


Fig.  13.  Jet  from  a  slot  of  uniform  velocity  and  adjacent  still 

width  a  fluid;  and  the  zone  of  establishment  of 

the  turbulent  jet  issuing  from  a  slot  of  definite  width  as  shown  in  fig.  13 
indicated  later,  he  considered  to  be  made  up  of  two  of  the 
vted  by  Tollmien  and  shown  in  fig.  11. 
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Albertson,  et  al.,  and  corresponds  essentially  to  the  second  aspect  of 
free  turbulence  treated  by  Tollmien— that  is,  the  efflux  of  a  jet  from  a 
two-dimensional  linear  slot. 

66.  In  this  part  (II),  To linden's  analysis  of  free  turbulence  at  a 
plane  of  velocity  discontinuity  and  Ms  analysis  of  the  lateral  diffusion 
of  turbulent  mixing  for  a  two-dimensional  jet  issuing  from  a  linear  slot 
will  be  applied  to  the  problem  of  determining  the  streamlines  and  com¬ 
puting  the  velocities  in  the  flow  field  of  a  stream  of  water  issuing  from 
a  channel  into  a  large  body  of  quiet  water  such  as  an  ocean. 

General  Considerations 


67.  Two-dimensional  turbulent- jet  theory  is  amply  treated  in  the 
literature,  and  for  the  pur-pose  of  this  report  there  is  no  need  to  repeat 
the  detailed  argument  here.  However,  as  will  be  seen,  certain  portions  of 
the  theory  are  necessary  for  proper  application  of  Tollmien1  s  computed  re¬ 
sults  to  the  problem  treated  herein,  and  for  this  reason  a  portion  of 
Tollmien' s  development,  with  only  minor  modifications,  will  be  sketched. 
Zone  of  flow  establishment 

68.  The  assumption  of  similarity  of  the  flow  characteristics  from 
section  to  section  of  the  jet  is  made  in  most  developments  of  jet  theory, 
and  has  been  amply  verified  by  experiment.  This  assumption  may  be  stated 
thus: 


(W) 

o 

where  u  is  the  horizontal  component  of  the  velocity,  U  is  the  velocity 

o 

in  the  undisturbed  stream  of  fig.  11,  b  is  the  width  of  the  diffusion 
zone,  and  0  is  a  constant  wMch  may  for  the  moment  he  considered 
arbitrary.  The  origin  of  the  coordinate  system  is  as  shown  in  fig.  11. 

69.  Tollmien  assumed  from  observation,  and  Liepman  and  Laufer  ^  P  ^ 
showed  analytically  that 

b  *  CSC  (49) 


where  a  is  a  constant.  The  assumption  of  similarity  of  flow  in  the  dif¬ 
fusion  area  thus  becomes 
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where  £  = 


r-'(') 

o 


(50) 


The  stream  function  is  determined  "by  using  equation  50: 


*  -  /udy  =  Uq  x  £F(n)  (51) 

where  F(tj)  =  /f (rj)drj  .  From  the  above  relation  the  vertical  component  of 
the  velocity  is  determined. 

v  - -  ~uo£F(n)  +  uo£nF'(Ti)  (52) 


Tollmien  assumed  the  Brandtl  expression  for  the  apparent  shear  stress, 
which  since  c*u/dy  is  positive  for  the  velocity  distributions  of  fig.  11  is 
given  by 


T 


(53) 


where  p  is  the  density  of  the  fluid,  and  -i  is  the  mixing  length. 
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70.  As  Goldstein  ’  *  7  has  pointed  out,  it  can  he  shown  from  the 

assumption  of  similarity  that  -t  must  he  proportional  to  h  ;  and  in  like 

X0  "o  Q 

manner  it  has  heen  noted  hy  both  Goldstein  and  Tollmien  ’  on  the 

basis  of  the  requirement  that  both  sides  of  the  equation  of  motion 


du 


+  v 


du 
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dr 

$y 


(5*0 


must  vary  in  the  same  manner  with  x  ,  that  -t  must  he  proportional  to 
x-~that  is, 

't  =  cx  (55) 

where  c  is  a  constant. 

71.  Substitution  of  equations  50,  52,  53#  and  55  into  equation  5^ 


yields 


2 

_F(ti)  =  F,n  (ti) 

r 


(56) 


From  equation  50,  t,  -  a&,  where  S  is  an  arbitrary  constant.  For  conveni¬ 
ence  (3  may  be  assigned  a  value  such  that 
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$3  «  2c2  (57) 

and  equation  56  becomes 

F(tj)  +  F'"  (n)  =  o  (58) 

which  is  the  differential  equation  derived  by  Tollmien. 

72.  The  boundary  conditions  follow  from  the  fact  that  at  the  upper 
boundary  of  the  diffusion  area  of  fig.  11,  u  *  U  and  du/dy  =  0,  and  at 
the  lower  boundary  u  =  0  and  bufby  =  0.  Letting  tj^  be  the  value  of  tj 
along  the  upper  boundary,  and  t]0  the  value  of  T|  along  the  lower  boundary, 
it  follows  that 

F  (tj  )  =  n 

F'(ni)  «  1 

F”^)  «  0  (59) 

F*(tj2)  =  0 
and  F,!(tj2)  =  0 

73.  Equation  58#  with  the  boundary  conditions  of  equation  59>  has 
been  solved  by  Tollmien,  and  the  results  are  shown  in  table  5. 

74.  In  Tollmien* s  solution  of  the  jet  boundary  problem  of  fig.  11, 
the  constant-velocity  stream,  as  is  evident  from  the  boundary  conditions  of 
equation  59,  is  assumed  to  be  unaffected  by  the  turbulent  mixing  process 
for  values  of  q  greater  than  tj^.  Hence  the  assumption  that  the  zone  of 
flow  establishment  of  fig.  13  consists  of  a  constant -velocity  inner  core 
with  a  diffusion  zone  on  either  side  corresponding  to  Tollmien  *  s  problem 
appears  to  he  warranted. 

Zone  of  established  flow 

75*  Tiie  problem  of  two-dimensional  efflux  from  a  linear  slot  was 
solved  by  Tollmien  in  a  manner  similar  to  his  solution  for  the  jet 
boundary.  As  before,  the  assumption  of  similarity  of  flow  from  section 
to  section  of  the  jet  is  made.  Equation  has  been  shown  by  IAepman  and 

7  p  12 

Lauf er  3  *  to  be  valid  for  the  zone  of  established  flow  as  well  as  for 
the  jet  boundary,  and  hence  the  condition  of  similarity  becomes 

=  f(tj)  TJ  =  (60) 

umax  * 

76.  Letting  U  =  U  for  x  *  x  ,  it  follows  directly  from  equation 
max  00 

60  and  the  momentum  equation. 


u  dy 


x  =  x 


that  IJ  cx  =  U  x  and  equation  60  then  becomes 
o  o  max 


if(,)  -y5**(ii) 


where  F’  (4)  *  f  (q)* 

77.  From  equation  62  the  stream  function  is  determined  to  be 


*=r  udy  *  U  J7  j;  £f(tj) 

w  O  o 


(63) 


and' 


•y  «« 


-  qF’(q)] 


(64) 


Referring  to  fig.  12,  it  is  apparent  for  the  upper  half  of  the  jet  where  q 
is  positive  that  du/dy  is  negative  and  hence: 


_  j>2  ,chi\ 

T  *  "P  ^ 


(65) 


78.  As  with  the  jet  boundary,  t  must  be  proportional  to  b  and 
x  ,  and  substituting  equations  55,  62,  64,  and  65  into  the  equation  of 

Q  2 

motion,  equation  54,  and  letting  as  before  =  2c  ,  we  have 


[f'(tj)J2  +  F(q)F"(q)  =  2F,,(q)Fm(q) 


id  direct  integration  with  respect  to  4  yields 
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F’(qMq)  =[F"(q)lJ 

M.  «» 


which  is  the  differential  equation  of  Tollmien. 

79.  The  boundary  conditions  follow  from  the  fact  that  at  q  =  0, 
u  *  U _ and  v  =  0,  and  at  the  boundary  for  q  =  q  ,  u  =  0.  From  which  it 


F(0)  =  0  and  F' (q_ )  *  0 


1< 


8l.  Since  the  appearance  of  Tollmien’s  paper  in  1926  there  have  been 
many  experimental  investigations  of  the  diffusion  of  submerged  two- 
linens ional  jets.  Without  exception,  insofar  as  the  writer  is  aware,  these 
investigations  have  been  made  with  air  jets  issuing  from  relatively  small 
slots.  Because  of  the  laws  of  similitude,  it  is  possible  to  apply  these 
experimental  results  with  confidence  to  water  jets  of  comparable  Reynolds 
numbers.  However,  owing  to  the  use  of  air  as  the  experimental  medium  and 
to  the  use  of  small  slots,  the  range  of  Reynolds  numbers  covered  by  the  ex¬ 
perimental  data  is  substantially  below  that  which  might  prevail  in  the 
present  problem  with  channel  widths  on  the  order  of  100  f fc  or  larger,  lor 
this  reason  it  became  desirable  to  establish,  from  previous  experimental 
work,  the  effect  of  Reynolds  number  on  the  constant  "c"  of  To  linden's 
analysis. 

Determination  of  Toll-mien's  "c” 

82.  gone  of  established  flow. 

,4  _1 _ _ _ _ _ In.  order  to  establish  the  relation  o 

V  Tollmien '  s  c  to  the  Reynolds  numbe 

\  for  the  zone  of  established  flow,  th 

M - - value  of  c  has  been  computed  from 

_ the  experimental  data  of  Albertson, 
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Dai,  Jensen,  and  Rouse,  Forthmann,' 
and  Reichardt.^  Hie  values  of  c 
thus  obtained  have  been  plotted  in 
fig.  l4  as  a  function  of  the  Reynolc 
number  R  =  U  a/v.  Inasmuch  as  the 

if*  * 


curve 
ical  j 


indicated  "by  the  Tollmien  analysis  for  velocities  near  the  center  line  of 
the  jet  and  lower  for  velocities  near  the  jet  boundaries. 

83.  Because  of  this  fact  the  problem  of  determining  Tollmen's  c 
from  experimental  data  involves  determining  the  best  fit  of  various  theo¬ 
retical  curves  to  an  experimentally  determined  curve  of  a  slightly  dif¬ 
ferent  shape.  In  order  to  avoid  the  indefiniteness  inherent  in  such  a 
process  it  appeared  desirable  to  define  the  experimentally  determined  c 
as  that  constant  which  caused  the  theoretically  determined  Tollmien' s  ve¬ 
locity  distribution  curve  to  coincide  with  the  experimentally  determined 

curve  at  u/u  =  0.5.  This  is  a  method  of  evaluating  constants  pertaining 
v  max,  ,,  ,  5,11 

to  jet  flow  that  has  "been  used  by  several  investigators. 

84.  In  addition,  Fortbmann,  in  determining  the  origin  of  the  ideal¬ 
ized  two-dimensional  jet  of  fig.  12  with  respect  to  the  two-dimensional 
slot  from  which  the  physical  jet  issues,  used  a  refinement  regarding  the 
location  of  the  origin  with  respect  to  the  slot  which  somewhat  complicates 
the  procedure,  and  which  was  not  demanded  by  the  accuracy  of  his  experi¬ 
mental  work.  In  order  to  simplify  the  procedure  of  computing  Tollmien  «s  c 
from  the  various  experimental  investigations,  the  origin  of  the  coordinates 
in  which  the  flow  distribution  is  expressed  was  taken  as  the  intersection 
of  the  center  line  of  the  jet  with  the  face  of  the  two-dimensional  slot. 
This  procedure  for  the  zone  of  flow  establishment  will  not  produce  ap¬ 
preciable  error  for  moderate  and  great  distances  from  the  slot.  This 
method  of  determining  Tollmien* a  c  has  the  advantage  of  being  simple  and 
definite.  In  fig.  15,  as  an  illustration,  Fortbmann* s  data  on  the  x~ 
component  of  the  velocity  has  been  plotted,  together  with  Tollmien* s  theo¬ 
retical  curve  for  c  «  0.0221*  (^2?  -  O.l).  It  will  he  observed  that  the 
experimental  theoretical  curves  coincide  at  u/u^  =  0.50*  Is  connec¬ 
tion  with  the  value  of  c  herein  determined  from  Fortbmann  *  s  data,  it  is 
proper  to  note  that  Forthmann,  on  the  basis  of  the  same  essperimental  data, 

_ value  almost  identical  with 
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Fig.  15.  Forthmann’s  data  on  air  Jets 


that  determined  by  Tollmien  for  the  free  Jet  boundary. 

85.  In  fig.  14  the  ■value  of  Tollmien1 3  c  determined  from  Forth- 
rnann's  data  for  the  zone  of  established  flow,  as  well  as  values  of  c  de¬ 
termined  in  a  similar  manner  from  the  data  of  Reichardt  and  Albertson  et  al., 
have  been  plotted  against  the  Reynolds  number  ^  .  Ua/v.  In  none  of  these 
experimental  investigations  was  the  air  temperature  given,  and  in  computing 
the  Reynolds  number  of  the  flow  the  air  temperature  was  assumed  to  be  70  F. 

B6»  From  fig.  14  it  is  evident  that  for  small  Reynolds  numbers, 
Tollmien's  c  decreases  rapidly  with  increasing  Reynolds  number,  and  as 
would  be  expected,  it  appears  that  the  value  of  c  approaches  an  asymp- 
totxc  value  as  continues  to  increase  to  the  poxnt  where  viscous  effects 
become  negligible.  The  data  presented  in  this  figure  are  obviously  some¬ 
what  fragmentary,  and  there  is  appreciable  scatter.  In  particular,  the 
curve  for  the  higher  values  of  the  Reynolds  number  is  defined  only  by  the 
experimental  data  of  Forthmann,  and  inasmuch  as  the  Reynolds  number  ob¬ 
taining  with  the  present  problem  of  flow  from  a  channel  into  an  ocean  or 


into  a  lagoon  may  be  on  the  order  of  1000  times  greater  than  those  shown 
in  fig.  14,  it  is  apparent,  because  of  gross  extrapolation  of  the  existing 
fragmentary  experimental,  data,  that  the  appropriate  value  of  c  for  the 
present  problem  is  somewhat  questionable.  Nevertheless,  from  the  data  of 
fig.  Ik  it  appears  that  the  value  of  c  approaches  asymptotically  a  value 
of  approximately  0.022  as  continues  to  increase  beyond  the  range  of  ap¬ 
preciable  viscous  effect,  and  it  is  proposed  to  use  this  value  of  c  in 
applying  To llmlen's  analysis  to  the  zone  of  established  flow. 

8?.  The  need  for  a  systematic  experimental  investigation  over  an  ex¬ 
tended  range  of  Reynolds  numbers  to  define  clearly  the  relation  between  NR 
and  c  or  some  other  comparable  parameter  of  the  jet  flow  such  as  is 
obvious. 

88.  Zone  of  establishment.  In  the  zone  of  establishment,  the  geom¬ 
etry  of  which  is  defined  in  figs.  11  and  13,  there  should  also  be,  at  least 
at  low  Reynolds  numbers,  a  region  where  viscous  effects  are  appreciable  and 
become  apparent.  In  this  connection  it  may  be  observed  that  there  is  no 
characteristic  constant  length  to  be  associated  with  the  geometry  of  the 
bounding  surfaces  of  the  jet,  and  therefore,  there  can  be  no  general 
Reynolds  number  defining  the  flow  characteristics,  which  remains  constant 
throughout  the  flow  pattern.  A  local  Reynolds  number  based  on  distance 
from  the  two-dimensional  slot,  width  of  the  diffusion  zone,  or  some  other 
length  parameter  of  the  flow  must  furnish  the  characteristic  length  of  the 
flow  phenomena  for  computation  of  the  Reynolds  number. 

89.  In  fig.  16  the  value 
c  stated  by  Tollmien^  and 

g 

Cordes  on  the  basis  of  their 
experimental  investigations, 
together  with  the  value  of  c 
computed  in  the  manner  described 
previously,  based  on  the  experi¬ 
mental  work  of  Idepraan  and 

7  2 

Laufer, ‘  Albertson  et  al. ,  and 

9 

Reicbardt,  have  been  plotted 
against  the  local  Reynolds 
number  U  x’/v,  where  x*  is 


Fig.  16.  Variation  of  Tollmien’s  c 
with  local  Reynolds  number  at  x'  =1  ft, 
zone  of  flow  establishment 
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the  distance  from  the  slot,  and  for  convenience  in  plotting  the  data  of 
fig.  16,  x"  has  been  taken  to  he  1  ft.  The  scatter  in  the  data  is  such 
that  a  precise  evaluation  of  c  is  impossible,  and  it  appears  reasonable 
to  assume  that  the  data  of  fig,  l6  is  in  the  region  where  the  Reynolds 
number  is  sufficiently  high  to  make  viscous  effects  negligible  and  hence 
to  make  c  Independent  of  Reynolds  number.  On  the  basis  of  the  data  of 
fig.  16  it  is  proposed  to  use  a  value  of  c  of  0.01?  in  the  application 
of  Tollmien1 s  results  to  the  zone  of  establishment. 

Determination  of  streamlines 

90.  i  establishment .  The  location  oi  the  ^}et  streamlxnes  in 
the  zone  of  establishment  may  be  determined  directly  from  Tollmien '  s  work. 
From  equations  52  and  58, 

* = uqx  F(ti)  (69) 

and  Tbllmien  has  computed  the  function  F(rj)  for  various  values  of  T|  . 
These  results  are  shown  in  table  5. 

91.  Measuring  x  in  terms  of  a  ,  the  width  of  the  slot,  we  may 
write  x  =  ka,  and  substituting  in  equation  69, 


Letting  Q  =  aUQ,  where  Q  is  the  volume  rate  of  discharge  through  the  slot 
per  unit  depth. 


(71) 


92*  From  equation  71  and  table  5?  k  =  x/a  may  be  computed  for  given 
¥’8.2x108  Of  and  i]  .  Using  the  value  of  given  in  the 

previous  section  we  may  determine,  from  the  assumed  tj  ,  the  value  of  y/a 
corresponding  to  the  x/a  =  k  and  t/Q  in  question.  This  has  been  done,  and 
the  results  have  been  plotted  in  fig.  17  to  form  the  streamlines  for  the 
zone  of  establishment.  In  this  connection  it  is  obvious  from  fig.  13  and 
the  constants  regarding  jet  diffusion  determined  by  Tollmien  (table  5)  that 
the  length  of  the  zone  of  flow  establishment  xq  may  be  determined  from  the 
relation 
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96.  Knowing  P(tj)  and  i]  from  Tollmien’s  computations  in  table  5  and 
using  the  value  of  c  =  0.022  or  0.0989  from  page  33,  the  values 

x/a  for  any  streamline  may  be  computed  from  equation  75*  Knowing  x/a  and 
rj  ,  the  corresponding  value  of  y/a  may  be  computed  Prom 


(76) 


97*  At  this  point  it  is  necessary  to  point  out  a  difficulty  which 
arises  when  Tollmien’s  analyses  for  the  jet  boundary  of  fig.  11  and  the 
two-dimensional  jet  flow  from  the  linear  slot  of  fig.  12  are  used  to  com¬ 
pute  the  characteristics  of  a  two-dimensional  jet  issuing  from  a  slot  of 
definite  width.  It  is  obvious  that  the  momentum  in  the  zone  of  flow 
establishment  must  be  identical  with  the  momentum  in  the  zone  of  estab¬ 
lished  flow.  Whatever  use  is  made  of  Tollmien’s  analyses  must  he  such  that 
this  elementary  momentum  consideration  is  observed.  In  fig.  3B  is  shown 
the  velocity  distribution  in  the  two  zones  of  the  jet  for  a  value  of 


Fig.  iB.  Velocity  distribution  at  x  =  x  =  6.12a 

o 
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X  .  X  =  6.12a,  the  nominal  limit  of  the  zone  of  establishment  as  given  in 
o 

the  preceding  section. 

98.  The  curve  for  the  zone  of  established,  flow  in  fig.  18  has  been 
computed  on  the  basis  that  U  =  U  at  x  =  6.12a,  and  both  velocity  dis- 
tribution  curves  have  been  computed  with  the  values  of  c  derived  in  par. 
82-87.  From  this  figure  it  is  apparent  that  the  momentum  consideration 
stated  above  will  not  hold  unless  an  adjustment  of  some  nature  is  made. 

This  adjustment  may  conveniently  be  made  by  considering  the  jet  to  consist 
of  three  zones  of  flow,  instead  of  the  two  heretofore  mentioned.  The  first 
zone  would  be  the  zone  of  flow  establishment  in  which  the  diffusion  of  the 
eddies  from  the  jet  boundary  steadily  expands  into  the  constant-velocity 
core  of  the  jet  until  the  core  is  annihilated.  The  second  zone  would  be 
the  zone  of  completely  established  flow,  and  the  third  zone  would  be  a 
transition  zone  between  the  first  and  second  in  which  neither  of  the  anal¬ 
yses  of  Tollmien  would  be  valid. 

99.  The  condition  that  the  momentum  in  the  zone  of  established  flow 
be  equivalent  to  that  in  the  zone  of  flow  establishment  requires  that 


(77) 


and  since 


=  f(n)  «  F*(n)  i)  =  ■  x  (?8> 

max  '\y  pc2 

we  have  by  equation  77 > 

Uo2a  -  2x  \^2?  F'(n)dn  -  xU^I  (79) 

where  _ _  /•»«> 

1  =  2  V2c2|  F*(t])dTj 
J0 

100.  Referring  to  the  idealized  jet  configuration  shown  in  fig.  13? 
in  which  the  zone  of  established  flow  begins  immediately  after  the  zone  of 
flow  establishment,  it  is  apparent  that  at  x  *  xQ,  -  Uq.  Using  this 

fact  in  equation  79  we  have 


101.  The  integral  I  may  he  evaluated  graphically  from  Tollmlen'f 
computations  and  is  found  to  he  0.1376.  Hence 


-f  ■  <dhr  -  T-27 


102-  This  is  the  value  of  xq  for  the  zone  of  established  flow  re¬ 
quired  for  constancy  of  momentum  in  a  two-dimensional  jet  Issuing  from  a 
slot  of  width  a  ,  whose  velocity  distribution  is  that  derived  by  Tollmien 
for  a  two-dimensional  jet  issuing  from  a  linear  slot. 

103.  The  value  of  x^/a  *=7.27  derived  above  for  the  zone  of  estab¬ 
lished  flow  is  to  be  compared  with  the  value  of  x  /a  =  6.12  derived  in  the 

o 

preceding  section  for  the  zone  of  establishment .  Since  in  the  derivation 

of  the  value  of  x  for  the  zone  of  established  flow  it  was  assumed  that 
o 

tT  „  =  U  at  x  -  x  =  7-2?&,  while  the  value  of  x  derived  for  the  zone  of 
max  OO  o 

flow  establishment  indicates  that  U  is  less  than  U  for  values  of  x 

max  o 

greater  than  6.12a,  it  is  apparent  that  the  velocity  distribution  developed 
by  Tollmien  demands,  for  constancy  of  momentum  from  section  to  section,  the 
existence  of  a  transition  zone  between  the  zone  of  flow  establishment  and 
the  zone  of  established  flow. 

104.  Substituting  the  value  of  xQ/a  =  7-27  obtained  above  in  equa¬ 
tion  75,  rearranging,  and  letting  Q  =  aU  ,  where  Q  is  the  volume  rate  of 

O 

flow  thorough  the  slot  per  unit  depth,  we  have 


I - i - - - -  -  F(q) 

■FTzj-nji 

and  using  the  value  of  \j2c^  =  0.0989  as  previously  determined  for  tl 


i  =  0.266  4x/&  F(tj) 


4l 


of  y/a  haw  been  computed  from  the  relation  y/a  «  x/a  ij  and  haw 

teen  plotted  in  fig.  17. 

106.  For  the  streamlines  f/Q  3  0  to  t/Q  =  0.5,  there  is  very  close 
agreement  at  the  section  x  =  6.12a  between  the  streamlines  conrputed  for 
the  zone  of  establishment  and  the  zone  of  established  flow.  The  divergence 
between  the  streamlines  computed  for  the  two  zones  is  more  marked  for  the 
streamlines  t/Q  =  0.6  and  0.7,  and  as  will  be  noted  from  fig.  17,  portions 
of  these  streamlines  consist  of  a  dashed  line  indicating  that  interpolation 
was  necessary, 

107.  With  regard  to  the  streamline  +/Q  =  0.7,  its  intersection  with 
the  jet  boundary  is  dictated  by  the  lateral  velocity  distribution  along  the 
jet  boundary.  From  table  5,  the  velocity  of  inflow  across  the  jet  boundary 
for  the  zone  of  establishment  is 

-  0.379  ^2?  =  0.0316  (84) 

o 

108.  In  like  manner  from  table  5  the  velocity  of  inflow  across  the 
jet  boundary  in  the  zone  of  established  flow  is 

jp  *  -0.498  Jxfi  *  -  — =r  (85) 

O  **  3C/ 

In  fig.  19,  equations  84  and  85  have  been  used  to  plot  the  solid  portion  of 


cyrvs  shown  .  The  doi/ii^d  portios  of  the  curve  is  &rbx1#rEry» 

The  actual  velocities  in  this  transition  zone  are  unknown,  and  the  straight 
dotted  line  in  this  region  has  been  drawn  merely  so  that  the  lateral  veloc¬ 
ity  distribution  of  fig.  19  is  consistent  with  the  streamlines  t/Q  =  0.8 
and  above— that  is,  so  that  the  relation 


* 

Q 


°-5  +  5  fy* 


(86) 


is  satisfied. 

109.  The  intersection  of  the  f/Q.  =  0.7  streamline  with  the  jet 
boundary  was  determined  in  accordance  with  the  lateral  velocity  distribu¬ 
tion  of  fig.  19* 

110.  as  the  velocity  is  everywhere  tangential  to  the 
streamlines,  the  data  of  fig.  1?  will  give  the  direction  of  flow  at  any 
point.  Also,  since  the  velocity  at  any  point  in  a  two-dimensional  field 
of  flow  is  inversely  proportional  to  the  spacing  of  the  streamlines,  the 
data  of  fig.  17  could  he  used  to  determine  the  magnitude  of  the  velocity 
vector  as  well  as  its  direction.  Nevertheless,  it  will  he  more  convenient 
to  determine  the  velocity  hy  other  means. 


Determination  of  veloc¬ 
ities  along  streamlines 


111.  Zone  of  flow  establishment.  Tollmien  has  computed  values  of 
v/u  ^2^  and  u/lJ  =  F*(*i)*  and  these  are  shown  in  table  5-  lotting 


'  o 

we  have 


112.  The  ratio  of  the  magnitude  of  the  velocity  vector  V  to  the 
velocity  TJ  through  the  slot  has  been  computed  from  equation  87  and 
table  5  for° various  values  of  i\  ,  and  the  results  have  been  plotted  in 


^3 


Fig,  20,  Velocity  distribution  in  the  zone  of  flow  establishment 

manner  have  been  plotted  against  s/a  in  fig.  21,  where  s  is  the  distance 
along  the  streamline,  measured  either  from  the  face  of  the  slot  or  from  the 
jet  boundary,  whichever  is  appropriate.  In  the  case  of  f/Q  greater  than  or 


greater  than  0.5  the  initial  point  was  the  intersection  of  the  streamline 
with  the  jet  boundary. 

114.  Zone  of  established  flow.  Eliminating  I  from  equations  79 
and  80 : 


_  2  2 
U  x  =  U  x 

O  O  HEX 


Substituting  the  value  of  U  y  thus  obtained  into  equation  78, 


tT"  «  F'(ti) 


and  the  function  F*(ij)  has  been  computed  by  Tollmlen  and  is  given  in  table 


5.  In  like  manner  the  term  sTx/x"  — — — - 
Jetting  0  W 


has  been  computed  by  Tollmien. 


!ufc  +  v 


we  have 


r 

o 


(tr  +  (if  '^) 


Letting  V  2c  »  O.0989  as  determined  previously,  the  velocity  distribution 
of  fig.  22  has  been  computed  from  equation  91  and  table  5.  From  the  data 
of  figs.  17  and  22,  the  velocities  along  the  streamlines  have  been  plotted 
against  s/a  in  fig.  23,  where  as  before  s  is  the  distance  along  the 
streamline  from  its  intersection  either  with  the  face  of  the  slot  or  with 
the  nominal  jet  boundary  shown  in  fig.  17- 

115,  In  fig.  24  the  data  of  figs.  21  and  23  for  the  velocities  along 
the  streamlines  in  the  two  zones  have  been  combined  to  give  the  velocities 


essentially  two -dimens ional .  Obviously  this  will  not  always  be  the  case  in 
practice.  Further,  the  analysis  of  Tollmien  assumes  that  Prandtl's  mixing 
length  is  constant  over  each  lateral  section.  In  the  case  of  an  actual 
ocean  and  channel  this  cannot  be  true,  since  the  mixing  length  along  the 
bottom  of  the  ocean  and  at  the  free  surface  must  be  zero. 

117.  Xn  a  channel  of  more  than  nominal  length  connecting  a  lagoon 
and  ocean,  the  velocity  will  not  he  uniform,  as  has  been  assumed  herein, 
but  will  have  the  typical  horizontal  and  vertical  velocity  curves  of  turbu¬ 
lent  channel  flow.  The  action  of  waves,  littoral  currents,  and  density 
currents,  if  present,  may  substantially  modify  the  jet  velocity  distribu¬ 
tion  shown.  Under  these  circumstances  the  solution  given  must  he  consid¬ 
ered  at  best  as  only  a  first  approximation  to  the  velocities  actually  pre¬ 
vailing.  Nevertheless,  despite  the  obvious  sources  of  error  mentioned,  it 
is  believed  that  the  application  of  two-dimensional  jet  theory  to  the 
problem  of  determining  the  velocity  distribution  in  the  stream  of  water 
issuing  from  a  channel  into  an  ocean  will  provide  the  basic  first  step  for 
the  investigation  of  the  methods  of  controlling  the  shoaling  in  the 
vicinity  of  such  channel  entrances. 

Status  of  Problem 

llS,  Additional  research,  both  experimental  and  analytical,  is  war¬ 
ranted  on  the  problem  discussed  in  Fart  II.*  In  particular  the  variation  of 

Tbllmien's  c  or  the  distance  x  with  Reynolds  number  in  the  region  of 

0 

established  flew  should  be  determined  by  systematic  experimentation.  In 
this  connection  it  appears  certain  that  for  Reynolds  numbers  character¬ 
istic  of  the  flow  prevailing  in  nature  for  the  problem  herein  under  consid- 


be  desirable. 

119.  Further,  it  would  appear  that  the  conventional  jet  phenomenon 
as  used  herein  constitutes  only  a  temporary  phase  of  the  flow  picture  pre¬ 
vailing  in  the  efflux  of  water  from  a  channel  into  the  tranquil  water  of  an 
ocean.  To  be  specific,  suppose  that  the  connecting  channel  of  fig.  10  were 
originally  absent,  and  that  such  a  channel  were  constructed.  Initially  the 
flow  from  the  channel  would  be  such  as  to  conform  at  least  approximately, 
it  is  believed,  to  the  streamlines  shown  in  fig.  1?  and  with  the  velocities 
along  the  streamlines  shown  in  fig.  2^.  However,  it  is  obvious  from  the 
streamlines  of  fig.  17  that  deceleration  occurs  and  will  cause  deposition 
in  the  jet  area  of  the  suspended  or  bed  load  originally  eroded  from  the 
channel  and  lagoon. 

120.  Further,  since  the  water  in  the  jet,  as  is  evident  from  fig. 

17,  is  being  decelerated  in  the  inner  core  of  the  jet,  while  along  the 
boundaries  acceleration  is  occurring,  it  might  be  expected  that  deposition 
in  the  region  of  the  boundaries  might  be  at  a  slower  rate  than  in  the  inner 
core  of  the  jet,  despite  the  fact  that  the  acceleration  of  the  water  in  the 
region  of  the  jet  boundaries  is  occurring  because  of  the  lateral  transport 
of  momentum  from  the  1.-nr»*»r  core  with  consequent  lateral  transport  of  sus¬ 
pended  load.  In  this  connection,  reference  to  fig.  6  in  Fart  I 
indicates,  for  the  reverse  of  the  tidal  cycle  herein  considered— that  is, 
flow  from  the  ocean  into  the  channel— that  the  acceleration  of  the  water 
as  it  approaches  the  channel  entrance  is  greater  in  the  region  of  the  jet 
boundaries  than  in  the  internal  core  of  the  jet.  Ibis  fact  would  also  in- 


be  the  pattern  of  deposition  of  suspended  load. 

122.  It  is  believed  that  an  investigation  of  tbe  variation  in  the 
net  rate  of  deposition  of  suspended  load  implied  by  tbe  velocity  d  Z3~  S  lj  I' 
tions  given  in  Parts  I  and  II  herein  would  be  a  logical  next  step  in  the 
development  of  the  general  problem  of  shoaling  near  tidal  channels  or 
entrances. 
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Table  1 


Computed  Values  of  x/b  and  y/b  for  Various  Values  of  f/Q  and  «0/Q 

Two-Dimensional  Orifice 
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Table  2 


Computed.  Values  of  s/b  and  q/v  for  Various  Values  of  t/Q  and  n0/ Q 

Two-Dimensional  Orifice 
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Computed  Values  of  x/b  and  y/b  for  Various  Values  of  A Q 

for  a  Projecting  Entrance 
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pendix  applies  the  general  solution  to  a  specific  case.  The  case  selected 
assumes  the  following; 

Width  of  channel*  1000  ft 

Average  flow  velocity  in  channel*  2.0  ft  per  sec 

Depth  of  channel*  not  pertinent.  Depth  over  entire  flow  area 
assumed  uniform 

2.  The  flow  net  is  drawn  up  using  table  2  and  fig.  5  of  the  main 
report  as  described  in  the  following  paragraphs. 

3.  First*  compute  the  free  streamline  velocity*  ¥  *  from  the  formula 
(see  page  11  of  main  report); 

v  =  v yo.611 

where  ¥ffl  is  the  mean  velocity  in  the  channel.  ¥  in  this  case  is  the  ve¬ 
locity  of  the  free  streamline  after  it  enters  the  channel;  it  is  also  the 
velocity  of  all  streamlines  at  an  infinite  distance  downstream  from  the 
entrance  assuming  nonviscous  flow.  In  the  case  at  hand; 

¥  *  V yo.611  *  2.0/0.611 

¥  =  3*27  ft  per  sec 

4.  Hext  enter  table  2  and  convert  the  q/v  column  into  q,  by  multi¬ 
plying  the  q/v  value  by  3.27.  The  symbol  q  represents  the  velocity  at 
the  indicated  location.  (This  conversion  is  made  in  table  A1  of  this 
appendix. ) 

5.  Using  fig.  5,  plot  those  computed  values  of  q.  at  the  proper  in¬ 
tersection  of  the  ¥/Q  and  X0/Q  coordinates  as  shown  in  fig.  5* 

6.  Convert  the  x  and  y  grid  in  fig.  5  to  a  dimensioned  grid  by 
multiplying  the  nondimensional  x  and  y  coordinates  as  shown  in  fig.  5 
by  one-half  the  channel  width  (in  this  case  by  one-half  of  1000  or  500  ft). 
The  results  of  this  conversion  sire  shown  in  fig.  A1  of  this  appendix. 

7.  The  resulting  diagram  is  a  plot  of  the  flow  velocity  and  dlrec- 
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found  near  the  outer  boundaries  of  the  channel  rather  than  in  the  center 
of  the  channel  in  this  type  of  solution. 


Fig.  Al.  Example  of  velocity  distribution  at  entrance  to  a  two-dimensional 
channel.  Assumed  width  of  channel  is  1000  ft  and  assumed  mean  velocity 
through  channel  is  £.0  ft/sec.  (Note:  Only  one-half  of  channel  and  flow 

pattern  is  shown.  } 
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APPENDIX  B:  APPLICATION  OF  AN  EXPANDING  FLOW  NET 

TO  A  SPECIFIC  CONDITION 

1.  Part  II  of  the  main  report  gives  a  general  solution  to  the 
problem  of  direction  and  velocity  of  flow  leaving  a  confined  channel.  The 
problem  is  essentially  that  of  a  two-dimensional  jet.  This  appendix  ap¬ 
plies  the  general  solution  to  a  specific  case.  The  case  assumes  the 
following: 

Width  of  jet  at  exit,  500  ft 

Average  flow  velocity  in  channel  before  exit,  3*0  ft  per  sec 

Depth  of  channel,  not  pertinent.  Depth  over  entire  area  of  flow 
assumed  uniform 

2.  The  flow  net  is  drawn  by  first  referring  to  fig.  IT  of  the  main 
report.  The  coordinates  in  fig.  1?  are  converted  to  dimensioned  quantities 
by  multiplying  the  x/a  and  y/a  coordinates  by  500  ft,  the  assumed  width 

of  the  jet  at  exit.  This  is  done  to  construct  the  coordinate  background 
of  fig.  Bl. 

3.  The  velocities  to  be  computed  will  fall  on  the  t/Q  streamlines 
shown  in  fig.  IT.  This  required  the  measurement  of  a  value,  s/a,  from 
fig.  IT;  s/a  is  the  distance  measured  along  these  streamline  curves  using 
the  dashed  line  as  the  zero  point  for  all  streamlines  greater  than  0.5,  and 
the  exit  point  of  the  jet  for  the  0.5  and  all  lesser  streamlines.  This 
distance  is  measured  and  its  value  is  determined  by  its  equivalent  length 
on  the  x/a  coordinate  scale.  Table  Bl  shows  the  compilation  of  the  s/a 
value  for  four  of  the  t/Q  streamlines  (1.4,  1.2,  1.0,  and  0.2). 

4.  Fig.  23,  which  shows  a  diagram  of  V/Uq  as  a  function  of  s/a 
values  along  the  streamlines,  ^/q,  is  then  utilized  (V  =  velocity  at 
selected  point  and  U  =  exit  velocity  of  jet).  For  each  s/a  value  on  the 
selected  */9  etreamliie  the  value  of  v/Bo  is  read  from  the  diagram  of 
fig.  23.  The  velocity  is  then  determined  hy  multiplying  the  V/U  value  hy 
3.0  ft  per  sec,  the  assumed  exit  velocity  of  the  jet.  This  process  is 
repeated  for  a  sufficient  number  of  the  t/Q  streamlines  to  give  the  de¬ 
sired  definition  of  the  flow  pattern.  The  results  of  these  steps  for  four 
selected  streamlines  are  shown  in  the  fourth,  fifth,  and  sixth  columns 

of  table  Bl. 


5.  These  steps  were  followed  in  constructing  the  velocity  flow  net 
shown  in  fig.  Bl.  The  applicable  values  of  the  x  sjei(3.  y  c  ooi*dxxi,sti  6  s 
of  fig.  Bl  are  shown  in  table  Bl,  and  are  obtained  by  multiplying  the  x/a 
and  y/a  values  in  the  table  by  500  ft,  the  assumed  exit  width  of  the  jet. 

6.  The  zones  of  flow  which  are  indicated  in  fig.  Bl  are  explained 
in  paragraphs  60  through  115  of  the  report.  In  effect,  these  zones  define 
the  transition  of  the  flow  from  a  streamlined,  unidirectional  jet  to  a 
dispersive,  turbulent  flow  expanding  into  the  larger  body  of  water.  The 
zone  labeled  "transition  zone"  in  fig.  Bl  is  not  specifically  defined  in 
the  text;  however,  it  represents  the  sector  between  the  end  of  the  "zone  of 
establishment"  at  x/a  =  6.12  and  the  beginning  of  the  "zone  of  established 
flow"  at  x/a  .  7.27. 
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